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Abstract. We find that universal three-body physics extends beyond the threshold 
regime to non-zero energies. For ultracold atomic gases with a negative two-body 
s-wave scattering length near a Feshbach resonance, we show the resonant peaks 
characteristic of Efimov physics persist in three-body recombination to higher collision 
energies. For this and other inelastic processes, we use the adiabatic hyperspherical 
representation to derive universal analytical expressions for their dependence on the 
scattering length, the collision energy, and — for narrow resonances — the effective 
range. These expressions are supported by full numerical solutions of the Schrodinger 
equation and display log-periodic dependence on energy characteristic of Efimov 
physics. This dependence is robust and might be used to experimentally observe 
several Efimov features. 
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1. Introduction 

Ultracold quantum degenerate gases near Feshbach resonances have been studied 
intensively in recent years [1]. The tunability of the interactions between the atoms, 
characterized by their s-wave scattering length a, makes ultracold quantum gases an 
ideal testbed to study a variety of few-body problems. Three-body systems in the 
strongly interacting regime |a| — oo especially interest people due to their importance 
in many areas in physics, including condensed matter, atomic, molecular and nuclear 
physics [2, 3]. 

Central among all ultracold three-body phenomena is the Efimov effect, which 
occurs when the scattering length a is much larger than the characteristic length scale 
tq of the atomic interaction [4, 5]. First predicted in the early 70's by the nuclear 
physicist Vitaly Efimov [6], this effect refers to the existence of an infinite number of 
weakly bound three-body states when \a\ — )■ oo. It has a profound impact on the three- 
body scattering observables and leads to universal three-body behavior irrespective of 
the short-range interactions between the particles [3]. In particular, when a < 0, the rate 
K3 for three-boson recombination, B + B + B ^ B2 + B, features a log-periodic series 
of resonance peaks on top of an overall scaling near zero energy. Explicitly [3, 4, 5], 
in atomic units, 

4590 sinh(2r/) 

i\3 = 7: 7) , (1) 

sin^[soln(|a|/ro)+$ + 1.53]+sinh ?7 m 

where sq ~ 1.00624 is a universal constant, $ is the three-body short-range phase [3, 7, 
8], ?7 is a parameter characterizing short-range inelastic transitions [3, 7], and m is the 
atomic mass. Physically, the peaks appear in when an Efimov state moves across the 
three-body break-up threshold and becomes bound [4, 9]. Another three-body inelastic 
process which is also important in ultracold experiments is the vibrational relaxation of 
a weakly-bound molecule: B2+B ^ B2 + B {a > 0). The rate for vibrational relaxation 
Vrei has the same functional form as Eq. (1), with the overall scaling changed to just 
a; the overall coefficient, to 20.03; and the constant in the phase, to 1.47. 

Near a narrow Feshbach resonance, however, the universal three-body physics is 
modified [10, 11, 12]. The dependence of the three-body physics on the resonance width 
in magnetic field can be parametrized with the two-body effective range res, which 
is inversely proportional to the resonance width [10]. It has recently been predicted 
that near a narrow resonance three-body inelastic processes dominated by short-range 
transitions are suppressed for bosons and enhanced for mixed spin fermions [12]. 

The Efimov and universal three-body physics introduced above are studied through 
the dependence of the three-body observables on a in the limit of zero scattering energy. 
Studies away from zero energy have mainly considered only the corrections to the zero- 
energy results [13, 14, 15, 16, 17, 18]. But, interestingly, it has recently been shown that 
the energy dependence also displays signatures of the Efimov effect for a > [19]. For 
energies beyond the threshold regime, however, the higher partial wave contributions 
to the inelastic rates are not suppressed by the Wigner threshold law [20] and tend 
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to mask these energy- dependent Efimov features. These contributions, combined with 
thermal averaging effects, make observing energy-dependent Efimov features in this 
case quite challenging. To overcome these difficulties, a scheme to collide Bose-Einstein 
condensates has therefore been proposed [19]. 

In this paper, we show that signatures of the Efimov effect also appear at non-zero 
collision energies in recombination for a < 0. In contrast to the finite-energy Efimov 
features identified previously [19], the Efimov signatures we show here have negligible 
higher partial wave contributions and are much less sensitive to thermal averaging. We 
also study the energy dependence of the vibrational relaxation rate when a > 0. Finally, 
we show that Efimov and universal features remain for large negative r^s, corresponding 
to a narrow Feshbach resonance. In most cases, we find analytic expressions for the rates 
and confirm them with numerical solutions of the Schrodinger equation. 



2. Method 



To study three-body universality at finite energies, we have numerically solved the 
three-body Schrodinger equation for identical bosons using the adiabatic hyperspherical 
representation [21]. In this representation, the Schrodinger equation reduces to a set of 
coupled equations for the hyperradial wavefunctions Fy{R), 

+ F, + Y,W.^'Fu' = EF,. (2) 

The hyperradius i? is a measure of the overall size of the system and ii = m/\/3 is the 
three-body reduced mass. The effective three-body potentials are 

W,,{R) = U,{R) - j^Qu.{R), (3) 

where Uy{R) are calculated from the adiabatic equation 

H^{R- ^)^,{R- fi) = U,{R)^,{R- fi). (4) 

The adiabatic Hamiltonian ifad includes the hyperangular kinetic energy and all 
interactions. The diagonal couplings Qyy{R) are determined from 



dR 



dR 



(5) 



where the double brackets indicate integration over the hyperangles Vt. The non- 
adiabatic couplings Wyi^i {v ^ v') govern the transitions between channels. The effective 
three-body potentials Wy^ are crucial for the study of three-body scattering, especially 
at low energies since they determine the universal behavior of the rates like Eq. (1) [7]. 

To solve Eq. (4) numerically, we must specify the interaction potential. Fortunately, 
when \a\ ^ ro, the universality of low-energy three-body physics allows us to model the 
interactions between the atoms near a Feshbach resonance by single channel interactions 
with large |a|, as the details of the atomic interactions affect three-body observables only 
through the short-range parameters $ and r]. This model is equivalent to a multi-channel 




Figure 1. The adiabatic hyperspherical potentials for a=-180 a.u. The deeply bound 
atom-molecule potentials are not shown, so all of the curves represent the three-body 
continuum, B+B+B, asymptotically. 



treatment as long as there is only one open channel in the two-body system. Therefore, 
when solving Eq. (4), we simply use a pair- wise sum of two-body potentials 

v{r) = — Dsech^(r/ro), (6) 

where r is the interparticle distance and D is a parameter that we use to control the 
scattering length as well as the number of two-body bound states. 

Our numerical approach has been described in detail in Refs. [21, 22, 23, 24]. Briefly, 
we first solve Eq. (4) by expanding $i/(-R; fl) on Wigner D functions for the Euler angle 
dependence of the overall rotations [21] and on a two-dimensional, direct product B- 
spline basis for the hyperangles representing relative motion in the body frame [22]. 
Inelastic rates are then calculated by solving Eq. (2) with the eigenchannel i?- matrix 
method [23, 24]. 

We use these numerical solutions to verify our analytical expressions for the 
universal behavior. These expressions are obtained by solving the hyperradial equation 
Eq. (2) using the universal results for the hyperradial potential W^,y{R) [25, 7]. Moreover, 
for recombination with a < and vibrational relaxation with a > 0, we need only 
know the behavior of the initial channel. In Fig. 1, we show W,^uiR) for a=-180 a.u. 
numerically calculated from Eq. (4). The lowest energy curve is the initial channel for 
recombination. It has a barrier at i? ~ |a| [26]. For smaller R, the potential has Efimov 
character; and for larger R, free particle character. For a < and \a\ ^ tq, we can 
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idealize the behavior in the incident channel for recombination as 

For R < To, W^i, is not universal and the solution from that region will simply be 
parametrized at R = tq with $ and rj. For vibrational relaxation when a > 0, the 
potential for the incident channel is only changed for i? ^ a to 

where E2b = is the binding energy of the weakly-bound molecule and / is the 

angular momentum between the atom and the molecule. 

For both processes, the coupling between the initial channel and the final deeply- 
bound atom-molecule channel is only significant for R < tq, and we can treat the 
channels as uncoupled for R> r^. The solutions are thus completely determined by the 
potential Eq. (7), and, by assuming each potential Eq. (7) holds over the entire domain, 
the hyperradial wave functions Fy{R) are [14] 

{Cosin($ + z?7) R = ro, 

CiR'/^Jis,{kR)-tiinSiNiso{kR)} ro < R < f3\a\, (9) 
C2R'/HJi^^+i/2{kR)-tRn6Ni^^+,/2{kR)} R > /3|a|, 

where J^ix) and N^{x) are Bessel functions of the first and second kind, respectively; and 
k = a/2/ZE with E the total energy relative to the three-body break-up threshold. Since 
the boundary between regions is not exactly at i? = \a\ and cannot be precisely defined, 
we include the constant /3 as a free parameter to be determined by matching to numerical 
solutions [12]. The effective angular momentum l^s is 3/2 for recombination [20] and 
for relaxation. The parameters $ and 77 are introduced in the same fashion as in 
Refs. [3, 7] to characterize the non-universal, short-range physics, including the inelastic 
transition. They will generally be taken to be free parameters and fit to the rates for a 
particular system. 

By matching the hyperradial wavefunction F^(R) at each boundary, the asymptotic 
phase shift 6 can be expressed in terms of P\a\ and the short-range parameters at any 
energy. To obtain the probability for inelastic transitions in such a single-channel model, 
we first calculate the probability for elastic scattering TZ which is given by the reflection 
coefficient: 



7^ 



2 

1 — itan5 



(10) 



1 + i tan S 

So long as 77 is not zero (which corresponds to no coupling to deeply-bound states), S 
will be complex and TZ will be less than unity. The overall probability for an inelastic 
transition is thus V = 1 — TZ. For three-body recombination of identical bosons, the 
recombination rate A'3 is related to V by [4] 
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Equation (1) can then be reproduced by deriving V for recombination near zero energy 

3. Three-body inelastic processes near broad Feshbach resonances 
3.1. Three-body recombination at finite energies 

In general, each region identified in Eq. (7) becomes important when the de Broglie 
wavelength becomes small enough to sample that length scale. Specifically, at energies 
above Eg = l/^r^, the details in the short-range region will start to play a role and 
the rates will not be universal. Below — that is, aX R > tq — the three-body 
potentials and couplings are all universal. We therefore expect universal behavior in the 
three-body scattering for energies below Eg. 

Because the Efimov region of the lowest potential in Fig. 1 and Eq. (7) has a barrier 
at i? I a I, it can support a three-body shape resonance when an Efimov state exists 
behind the barrier. When the energy of the Efimov state is above the three-body break- 
up threshold E=0, increasing the collision energy across the resonance energy will lead 
to a resonant peak in K3 [4, 9]. This can be understood from the fact that near the 
resonant energy, the amplitude of the hyperradial wavefunction behind the barrier in 
W^i, is greatly enhanced. The three-body system thus has much larger probability in 
the small-i? region where the coupling between the incident channel and the deeply- 
bound atom-molecule channel peaks. However, because the height of the barrier is 
~ E'b = the resonant behavior in can only be observed for energies below E]^. 

And, below this energy, the barrier can only confine a resonant Efimov state with size 
on the order of \a\, implying that at most one Efimov resonance can be observed as a 
function of energy for a fixed a. 

In Fig. 2, we show the evolution of the first Efimov resonance peak from finite energy 
to zero energy as \a\ increases (see also Ref. [9]). These recombination probabilities 
were calculated numerically. Interestingly, when a peak moves below zero energy, an 
oscillatory structure with one full period appears in V ioi E > E\^ as shown in Fig. 2(b). 

To understand the oscillations revealed in these numerical results, we employ the 
analytical model introduced in Sec. 2. Matching the hyperradial wavefunctions at finite 
energies gives 



n 



{J'2-iK)[J,,,-tim5iNi,,] \f{J^ + iN^)[Jiso~tan5iNis, 



iJ2-^N,)[Ji^-t^n6,NiJ y\iJ, + tN,)[Ji^-te.n6^Ni^' ^> ' ^^^^ 

The primes denote derivatives, and the Bessel functions should be evaluated at i? = P\a\. 
When the scattering energy satisfies /c ^ l/rg, tan 61 is given by 

sin((/) -I- 2/7 — i7rso/2) , , /, x 

tan^i^ T^Y' tan0 = tan <l>-soln A:ro +¥'o + vr/4 , 13 

cos(0 + 2?7 + 27rso/2) 
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Figure 2. The three-body recombination probability 'P{a < 0) for three identical 
atoms with the mass of Cs when a resonant Efimov trimer state is above the three- 
body break-up threshold [27]. (a) As \a\ increases, the resonant peak moves towards the 
break-up threshold, (b) Increasing \a\ further after the Efimov state becomes bound, 
the resonant peak disappears, leaving an oscillatory structure behind. 



with 

_ Re[r(zgo)] - Im[r(zso)] 

'""^°-Re[^(^.o)] + Im[^(z.o)]• ^''^ 
For identical bosons, v9o=— 0.154187r. Equations (12) and (13) give a complete universal 
expression for in terms of the scattering length and energy. This analytic expression 
is evaluated and plotted in Fig. 3. In (a), the evolution of the resonant peaks with 
energy can be clearly seen, but the oscillatory structure — which should appear at 
fixed, large | d I cis Sb function of — is essentially invisible due to its relatively small 
modulation compared to the resonant peaks. Figure 3(b) replots the data to better 
show the energy- dependent oscillations. Note that the fine structure parallel to the E^, 
line is an artifact of the discontinuity of our idealized potential in Eq. (7) at R = (3\a\. 
The connection between the Efimov resonances and the energy-dependent oscillations 
is evident in Fig. 3(b). The peak of an oscillation appears approximately at the energy 
where the zero-energy position of each Efimov resonance intersects the E'b line. The rate 
oscillates through a full period when an Efimov resonance moves below the three-body 
breakup threshold. Therefore, the number of bound Efimov states can be directly read 
from the number of full oscillations. 

Remarkably, when 77 ^ 1 (corresponding to a small probability for an inelastic 
transition at small R as is expected), Eqs. (12) and (13) reduce to a simple form in the 
energy range E\, <^ E <^ {kr^ <^\): 

P ^ 2sinh(7rgo) sinh(2?7) ^^^^ 

cosh(7rso) + sin[— 2so ln(A;ro) -I- 2$ — 2(y9o] ' 
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Figure 3. The universal dependence of A'3 on scattering length and energy for 
identical bosons with the mass of ^Li with J—0 as evaluated from Eqs. (11) and (12). 
(a) The rate scaled as K-^E^ /a^ to show better the resonant peaks, (b) The same data 
scaled as K^E^ to emphasize the energy-dependent oscillation described by Eq. (15) 
which can be seen in the upper right half of the plot. Note that the fine structure 
parallel to the E^ line is an artifact of the discontinuity of our idealized potential in 
Eq. (7) at i? = /3\a\. In both plots, <I>=7r/2, 77=0.02 and ro=15 a.u.. 

It should be noted that for E > the parameters $ and 77 will depend on energy 
in general. However, for energies E <^ the change in the shape of the short-range 
hyperradial wavefunction is small since E^ is still small compared to the depth of the 
potential, therefore both $ and 77 can be treated as energy-independent quantities. 

The log-periodicity of the energy-dependent oscillations in Eq. (15) with period 
g27r/so provides the clear connection with Efimov physics. Like the energy-dependent 
oscillations observed for a > [19], the overall phase of the oscillations does not depend 
on a as is clearly shown in Fig. 3(b). Therefore, as \a\ increases, the lower limit of 
validity E^, = 1/yua^ of Eq. (15) decreases and more oscillations will appear towards 
the zero-energy limit. At the same time, those that have already appeared at higher 
energies remain unchanged. 

These oscillations might thus provide one experimentally accessible means of 
measuring multiple Efimov features. In the experiments observing the Efimov effect via 
zero-energy three-body recombination [28, 29, 30, 31, 32, 33, 34, 35], the temperature 
must be kept below as |a| is increased [13]. This restriction makes the observation of 
multiple Efimov features quite challenging. However, if the Efimov features are measured 
via the energy-dependent oscillations, the experimental requirements might be less 
demanding. Furthermore, the a- independence of Eq. (15) allows the energy-dependent 
oscillations to be measured at different fixed a, provided that E > E^^ (/c|a| > 1). 
Because the modulation of the oscillations is small compared to the resonant peaks, 
though, the scattering length where the measurements are carried out should be chosen 
to be far away from resonant peaks. 
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Figure 4. (a) The three-body recombination probability V{E) for three identical Cs 
atoms, two identical Cs atoms with one Cs atom in a different spin state (Cs') and 
two identical Cs atoms with one Li atom. The symbols are numerical results, and the 
black solid lines are the analytical result from Eq. (15). The analytical curves are only 
shown over their range of validity, < E < E^, but E^ is off the scale. For CsCsCs, 
a = 2 X 10^ a.u.; for CsCsCs', a = 6 x 10^ a.u.; and for CsCsLi, a = 10^ a.u. (b) 
The recombination probability after thermal averaging. To make clear the comparison 
of the relative modulation in the oscillations from the three different systems, V for 
CsCsCs' and CsCsLi has been multiplied by 5 and 1.9, respectively. 



Our results apply equally well to three-body systems BBX with two identical 

bosons, assuming a ^ tq is the B+X scattering length. All that changes is the value 

of So [36]. The density of Efimov states in such systems grows with the mass ratio 

7 = niB/mx since sq grows with 7, decreasing the spacing e^'^/'^" of energy-dependent 

Efimov features [36, 25]. As can be seen from Eq. (15), however, the relative size of the 

modulation of V is 

, _ cosh(7rso) , , 

AP = 2 . J '\ (16) 
smh [tisq) 

which increases with decreasing sq. The two effects thus conflict: to see more Efimov 
features requires a large sq (7), but seeing them clearly requires a small sq (7). 
Figure 4(a) shows V oc E^K^ [see Eq. (11)] for model systems chosen to have mass 
ratios matching representative atomic systems. For three identical Cs atoms, Eq. (16) 
predicts AV to be 0.170; for two identical Cs atoms and one Cs atom in a different spin 
state, the ratio is 1.37; and for two identical Cs atoms with one Li atom, the ratio is only 
0.00789. The predicted spacing between features is 515, 3.9x10^, and 23.8, respectively. 
The analytic expression, Eq. (15), was fit to the numerical data by adjusting only $ 
and ?7 and is also plotted in Fig. 4. It agrees quite well with the numerical results. 

A practical balance of contrast versus spacing can likely be found by first 
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determining the smallest contrast experimentally observable. This choice makes sq as 
large as experimentally tolerable and thus gives the smallest spacing between features. 
For instance, if a relative modulation of A'P=10% can be observed, then from Eq. (16) 
we find sq = 1.1748. The corresponding mass ratio is 7 ~ 7.1, and the spacing between 
Efimov features is 210 — a considerable improvement over three identical bosons. 

To determine whether these features are experimentally observable even in the 
presence of a distribution of collision energies, we must thermally average the energy- 
dependent rates. Assuming a Boltzmann distribution [37, 13], we find the results shown 
in Fig. 4(b). For < E < Es, the oscillatory structure in V generally becomes less 
clear after thermal averaging. For BBX three-body systems where Sq is small, such as 
CsCsCs', the oscillations are well preserved after thermal averaging, with the trade-off 
that the oscillation periods are huge. 

3.2. Vibrational relaxation at finite energies 

Although the relaxation rates K-ci for three identical bosons with a > resemble the 
a < recombination rates in the zero-energy limit [3, 7], they differ dramatically at 
higher energies. Note that zero scattering energy now refers to the atom-molecule 
breakup threshold rather than the three-body breakup threshold as for recombination. 
The usual two-body Wigner threshold law gives a constant rate for atom-molecule 
collisions in the threshold regime E < i?th where generally i?th ~ -E'b- An exception 
arises when the zero energy value of Kei is near a resonant peak which occurs when the 
atom-molecule scattering length Oam is much larger than Tq. In this case E^ii ~ 1/ /iamOam; 
where /lam is the atom-molecule reduced mass. For i^th ^ E < Eg, K-ci scales like 

V^ci = Ck~\ (17) 

The constant C depends on the short-range physics and is related to the short-range 
parameter 77. But different from three-body recombination, the relaxation rates does 
not show any oscillations in our numerical calculations. 

4. Three-body inelastic processes near narrow Feshbach resonances 

Following Ref. [12] , we model a three-body system near a narrow Feshbach resonance 
by using pairwise two-body potentials that support an s-wave shape resonance. We 
accomplish this by adding a barrier to our model two-body potentials, ensuring that the 
resonance has a large, negative effective range Vcs [10, 11]. For such reg, the potential 
Wi,i,{R) in the region ro <^ -R ^ |rcff| is replaced by a weak, non- universal Coulomb 
potential [12]. For identical bosons, when |a| ^ |rcfr| ^ vq the new length scale produces 
a l/|refj| suppression in the zero-energy inelastic rates that lead to a deeply bound two- 
body state [12]. 

This new length scale in W^y{R) introduces another energy scale for the universal 
behavior when the wavelength is short enough to sample the three-body potentials and 
couplings in ro ^ -R ^ I'^cffl- The new energy scale is Ees = l//i|a'"eff|^ and defines a 
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new energy region: E^s < E < Eg. The parameter a ~ 0.28 plays a role similar to P 
and was determined in Ref. [12] by fitting the universal formula for near a narrow 
Feshbach resonance to numerical results near zero-energy. Although it was determined 
at zero energy, it should apply at non-zero energies as well. Where we treated these 
processes at zero energy in Ref. [12], we will now find the non-zero energy expressions 
for their rates. 



4.1. Three-body recombination at finite energies 

We calculate for identical bosons for a < numerically up to the short-range 
energy Eg. As shown in Fig. 5(a), the recombination probability V at E ^ Eg 
has only a relatively weak dependence on r^s and a, since in our model the short- 
range behavior does not have strong dependence on r^s or a. From our numerical 
calculations, we deduce that the energy dependence changes to P oc /c in the energy 
range Ecs < E < Eg. Therefore as the energy gets smaller, the recombination 
probability decreases monotonically until E = E^s, where oscillatory behavior similar 
to Eq. (15) takes over. Compared to a broad resonance (rcff=— 20 a.u. in the figure), 
V is suppressed by a factor of l/|refr| around E ^ E^s before it connects to the 
oscillatory behavior at lower energy. Considering that K3 at a fixed a has the same 
energy dependence for both broad and narrow Feshbach resonances for E < Ef,g, the 
recombination rates at lower energies are all suppressed by l/|rcfr|. 

To understand this observed l/|rcff| suppression, we follow the analytical procedure 
introduced in Ref. [12]. The Coulomb-like potential in the region |rcfr| ^ i? ^ |a| takes 
the form 

2/i|reff|i? 

with Co a non-universal constant on the order of unity that can be positive or 
negative [12]. We found in [12] that simply setting W^i, to zero gave analytic results 
consistent with our numerical results and is justified by the fact that W^i, is not universal, 
changing from attractive to repulsive, while the numerical rates are universal. To derive 
analytic expressions for the rates that connect most closely with those for a broad 
resonance, we write the wave function in the region tq ^ -R ^ |rcff| as 

F^ = Csm{kR + Ses). (19) 

For E <^ Eg, we can use a low-energy expansion for S^s, 

Sos ~ -Ak, (20) 

where we have introduced the complex three-body short-range scattering length A to 
parametrize the physics at R < tq [12]. The recombination probability V can thus be 
derived in the same fashion as the broad Feshbach resonance case. In the energy range 
E\j < E < Ecs, this gives the same results as Eq. (15), with tq replaced by a\rca\ and 
the parameters $ and rj connected to the real and imaginary parts of A: 

ImA 



^ ^ a — ReA/\res\ . , 

tan$ = 2so — , ^ . r, smh?7 

a + Re/l/|refr| 



sin 2$. (21) 




Figure 5. (a) The thrcc-body recombination probability 7^(0 < 0) for identical bosons 
near narrow Fcshbach resonances with roff=— 20, —200, and —2000 a.u.. For all cases, 
a=— lO"" a.u., and tq ~ 50 a.u.. The symbols are the numerical results, and the black 
solid lines are from Eq (15) and (21). The analytical curves are only shown over their 
range of validity, < E < E^s- Equation (22) is also shown for reS = —2000 a.u. 
with ImA fitted to be 1.1 x 10~* a.u. The parameters Re^ are directly calculated from 
the short-range three-body potentials under the single-channel approximation, giving 
Rcyl=50, 57, and 180 a.u. for r^ff = -20, -200, and -2000 a.u., respectively. The 
parameter a is then obtained by fitting to be 0.25, consistent with Ref [12]. (b) The 
recombination probability from numerical calculations after thermal averaging. 

These analytical results are shown together with the numerical data in Fig. 5(a). The 
deviation between the numerical and analytical results for rcff=— 20 and —200 a.u. is 
because the universal requirement |rcfj| ^ tq is not well satisfied. 

For Ees < E < Eg, the recombination probability V has a simple scaling behavior 
independent of the scattering length: 

V = 4:k\lmA\, (22) 

in agreement with the numerical results. 

To check the robustness of these oscillations for a thermal distribution, we show 
the thermally averaged V with large |rcfr| in Fig. 5(b). It can be seen that although 
the modulation of the oscillations at low temperature is reduced, the V oc k scaling for 
Ecs < E < Eg is essentially unchanged. 

4-2. Vibrational relaxation at finite energies 

The behavior of vibrational relaxation for identical bosons with a > changes in the 
same energy range as for three-body recombination: E^^ < E < Eg. Figure 6 shows 
that as |refr| increases, instead of the scaling seen in Eq. (17) for broad Feshbach 
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Figure 6. The three-body relaxation rate Vrd for identical bosons near narrow 
Feshbach resonances with a— 10^ a.u. and reff=-200, -1000. and -2000 a.u. The 
rates are shown for collision energies above the three-body break-up threshold, where 
no resonant peaks are present so that the energy dependence is more easily seen. The 
symbols are the numerical results, and the back solid lines are the thermally averaged 
results. 



resonances, a plateau region appears in Vre\, extending from Eg down to E^s. This 
change also connects to the l/|refr| suppression of the rates near zero energy and can be 
used for observing the effective-range effect near narrow Feshbach resonances. Further, 
this plateau behavior can be derived in the same manner as Eq. (22). 

Since the value of the plateau has only a weak dependence on rgg and extends down 
to Ecs, the larger |rcff | is, the broader the plateau will be. And, because Vrei for different 
Tefj has similar values for E ^ Eg, the fact that the plateau extends down to E^s will 
make K-ei smaller than for a broad resonance by a factor of ro/|rcff| for all energies 
E < Ecs- In other words, we expect C in Eq. (17) to be proportional to l/|rcff|. The 
experimental observation of a plateau region in the energy dependence of Vj-d will show 
the universal physics pertaining to a narrow Feshbach resonance. Moreover, the energy 
where V^ei changes from plateau to scaling can roughly give the value of |reff|. 

In Fig. 6 we also show the thermally averaged relaxation rates. It can be seen that 
the plateaus in E^^g < E < E^ are not dramatically changed by a thermal distribution. 

5. Higher partial wave contributions 

To determine whether the energy-dependent features we have described will be 
observable, the contributions of higher partial waves must be considered as they can 
easily obscure features away from threshold. Key to this consideration is the fact that 
three-body recombination (a < 0) and relaxation (a > 0) are dominated by transitions 
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at R < tq. For the three-body systems where the Efimov effect does not occur for J > 
— for instance, identical bosons and BBX systems with 7 < 38.61 [36] — the effective 
hyperspherical potentials are all repulsive for R> tq, leading to substantial suppression 
of the inelastic rates when E -C Eg. Following the same method used for deriving the 
energy dependence above, we find that the recombination probabilities for J > scale 
like (/cro)^^°, where po is a universal constant that increases with J [36]. For identical 
bosons, Pq is always larger than 2. For E <^ E^ [krQ <^ 1), V{J > 0) are all negligible. 
Moreover, we have confirmed this conclusion for a few cases numerically. For vibrational 
relaxation, there is similar suppression for J > contributions. 

6. Summary 

We have studied inelastic three-body collisions and have found that they behave 
universally for collision energies well out of the ultracold regime — up to nearly a 
milli-Kelvin in the examples shown. Moreover, we have identified much of the universal 
energy dependence as a manifestation of the Efimov effect, opening up new dimensions 
on Efimov physics. Significantly, we have been able to derive analytic expressions, which 
were motivated and verified by numerical calculations, for these rates in terms of a few 
parameters. And, we have done so for both broad and narrow resonances. Our study 
thus brings universal three-body physics out of the ultracold regime, suggesting even 
richer universal physics can be discovered both theoretically and experimentally. 
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